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Abstract —In the recent past, new enhancements based on the 
well established Aloha technique (CRDSA, CRDSA++, IRSA) 
have demonstrated the capability to reach higher throughput 
than traditional SA, in bursty traffic conditions and without any 
need of coordination among terminals. In this paper, retransmis¬ 
sions and related stability for these new techniques are discussed. 
A model is also formulated in order to provide a basis for the 
analysis of the stability and the performance both for finite and 
infinite users population. This model can be used as a framework 
for the design of such a communication system. 

I. Introduction 

Since its birth more than 40 years ago, ALOHA ||T| Q 
has established as one of the most well-known protocols 
for packet-switching wireless networks. Since then, several 
papers have been published on this random access technique. 
Some papers are concerned with the issue of improving the 
throughput (e.g. providing slots for synchronized transmission 
from users or providing diversity by sending the same 
packet more than once 0) while other papers concentrate 
on the study of the stability 0 0 and on the retransmission 
policies to ensure that the communication is stable ||7|. In fact, 
when retransmissions are considered, a feedback is created 
between the transmitter and the receiver. Therefore, depending 
on the arrival rate of packets to be transmitted and on the 
Packet Loss Ratio, the possibility of an overloaded or saturated 
channel is present. Eor this reason, studies on stability of such 
a communication scenario and on possible policies to have 
a channel working in optimal conditions have been of big 
interest. Recently, a new technique named Contention Reso¬ 
lution Diversity Slotted Aloha (CRDSA) has been introduced 
0 - This technique allows to reach values of throughput up 
to 0.55|^ using an approach similar to Diversity Slotted 
Aloha 0 (i.e. every packet is sent twice in different slots) and 
adding a Successive Interference Cancellation (SIC) process 
at the receiver in order to attempt restoring collided packets 
(details of how this is done will be illustrated in Section II). 
Afterwards, the same concept has been extended to more than 
two instance^per packet (CRDS Ah-h-) 0 and also to irregular 
number of packet repetitions (IRSA) As a consequence 
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* The term instances is used towards the paper with the meaning of ’’total 
number I of transmissions for the same packet in different slots”. 


of these new enhancements, also the study of the stability of 
these new techniques is of big interest. Therefore the objective 
of the study shown in this paper is the development of a 
model able to evaluate the stability for these new techniques 
when retransmissions are considered. To do so, we refer to 
0 in which a performance evaluation of Slotted Aloha was 
presented. In particular, the mentioned paper analyzes system 
stability behaviour through the development of analytic models 
based on the so called equilibrium contour, that represents 
those points for which the expected number of successful 
packets is equal to the average number of newly generated 
packets, so that the overall communication is stable (i.e. the 
expected number of packets sent is the same at any time). In 
the followings, we adapt this model to our case to provide a 
useful tool for the study of the stability in CRDSA, CRDSA-h- 
and IRSA. Moreover, a framework for the optimization will 
be introduced. The starting point and the aim of this paper 
are similar to those in m- However, while the graphical 
representation used by Kissling represents the drift of the 
communication as the one in 0, our representation is based 
on the one in 0. The details for this different representation 
choice will be clear throughout the paper, but the main point 
is basically the possibility to graphically distinguish a curve 
influenced by the retransmission probability and a load line 
determined by the number of users and their probability of 
transmission so that these dependencies are separated and 
it is easier to understand how the communication behaviour 
changes when changing those values. Moreover, in GD the 
goal of showing agreement between simulations and expected 
analytical drift is reached by averaging simulations results 
over a consistent number of trials for any possible initial 
state (in terms of number of backlogged users). In this paper 
instead, the result of single simulations with various settings 
are presented, aiming at showing the outcomes in terms of 
backlogged users, throughput and so on, frame after frame 
and starting from the initial state with no backlogged users. 
This allows to empirically understand how and how much the 
communication values move around the expected behaviour. 
Einally, an average packet delay analysis is also introduced. 

H. System Overview and Problem Statement 

Consider a multi-access channel populated by a total num¬ 
ber of users M (finite or infinite). Users are synchronized so 
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unsuccessful packet in one of the successive frames with 
probability pr. In the remainder of this paper, we analyze 
the conditions that ensure stability of the overall transmission 
when using such a policy and we give a simple yet effective 
tool to design such a communication scenario in relation to 
the expected throughput and packet delay distribution gathered 
from the model. 


Fig. 1: Example of frame at the receiver for CRDSA (2 instances 
per packet). Plain slots indicate that a transmission occurred for that 
user in that slot. 


that the channel is divided into slots and Nf consecutive slots 
are grouped to constitute a so called frame. The probability 
that, at the beginning of a frame, an idle user has a packet 
to transmit is pq. When a frame starts, users having a packet 
to transmit place I instances of their packet over the Nf slots 
of the frame. The number of instances I can be either the 
same for each packet (regular burst degree distributio^^ |[^ 
g or not (irregular burst degree distribution) fT^. Packet 
instances are nothing else than redundant copies as in Q 
except for the fact that each instance contains a pointer to 
the location of the other instances. These pointers are used in 
order to attempt restoring collided packets at the receiver by 
means of Successive Interference Cancellation (SIC). Consider 
Figure representing a frame at the receiver for the case of 
2 instances per packet (CRDSA). Throughout the paper, we 
assume perfect interference cancellation and channel estima¬ 
tion, which means that the only cause of disturbance for the 
correct reception of packet’s instances is interference among 
them. Moreover, FEC and possible power unbalance are not 
considered. Each slot can be in one of three states: 

• no packet’s instances have been placed in a given slot, 
thus the slot is idle; 

• only 1 packet’s instance has been placed in a given slot, 
thus the packet is correctly decoded; 

• more than 1 packet’s instance has been placed in a given 
slot, thus resulting in interference of all involved packets. 

If at least one instance of a certain packet has been correctly 
received (see User 4), the contribution of the other instances 
of the same packet can be removed from the other slots. 
This process might allow to restore the content of packets 
in slots where intereference occurred (see User 2) and in an 
iterative manner other packets may be correctly decoded, up to 
a point in which no more packets can be restored or until the 
maximum number of iterations for the SIC process is reached. 

Depending on the design choice, the packets that have 
not been decoded at the end of the SIC process are either 
discarded or a retransmission process is accomplished. In the 
latter case, a feedback is needed in order to inform users about 
the eventual failure of their transmission so that a certain 
retransmission policy can be applied. In the retransmission 
policy considered in this paper, each user retransmits its 

^The burst degree distribution is defined as the probability distribution to 
have a certain number of instances for a certain packet. 


III. Stability Model 


Consider the aforementioned random access communication 
system. Each user can be in one of two states: Thinking (T) or 
Backlogged (B). Users in T state are idle users that generate 
a packet in a frame interval with probability poi if they do, 
no other packets are generated until successful transmission 
for that packet has been acknowledged. Users in B state are 
users that failed transmitting their packet, therefore they are 
waiting either to retransmit their unsuccessfully transmitted 
packet with probability p^ at the beginning of each frame (thus 
geometrically distributed) or to receive a feedback about the 
outcome of their retransmission. Moreover, we assume that 
users are acknowledged about the success of their transmission 
at the end of the frame. 

Let’s define 

• Ng : backlogged packets at the end of frame j 


r<3 _ 

'^B — Nf 

to users in B 


— : expected channel load of frame j due 
state 


• Gip : expected channel load of frame j due to users in T 
state 


• Gjj^ = Gip + : expected total channel load of frame 

3 

• PLR3 {G^pj^, Nf,d, I max) '■ expected average packet loss 
ratio of frame j, with dependence on the expected total 
channel load Gin, the frame size Nf, the burst degree 
distribution d and the maximum number of iterations for 
the SIC process Imax 

. G^oup = G\pp{l - PLR{G^pjp,Nf,d,Imax)) : part of 
load successfully transmitted in frame j, i.e. throughput. 


Our aim is to find the equilibrium contour in the (Nb,Gt) 
plane defined as the locus of points for which at any frame, 
the expected channel load due to users in T state is equal to 
the expected throughput. Thus 

Gt = Gout = GiNi^ — PLR{Gin, Nf, d, Imax)) (1) 

where the frame number j has been omitted, since in equilib¬ 
rium state this condition must hold for any frame. Moreover, 
being in an equilibrium point implies that also the expected 
number of backlogged users remains the same frame after 
frame. Therefore 


NB=NB{l-Pr)+GlNNfPLR{GlN,Nf,d,Imax) (2) 

from which 

_ GinPLR{Gin, Nf, d, Imax)Nf 
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Fig. 2: Equilibrium contours for CRDSA with Nf = 100 slots, 
T =20 

-^max 


Equations Q and 0 completely describe the equilibrium 
contour. In fact, once the configuration parameters Nf, d, 
Imax of the system are set and the retransmission probability 
is chosen, the curves are completely described plotting Gt 
and Nb for different values of apart from the M 

and po considered. Figure displays examples of equilibrium 
contour for various values in the CRDSA case. As we can 
see, when the probability of retransmission on a given frame 
decreases, the equilibrium contour moves upwards so that the 
the maximum expected throughput is obtained for a bigger 
mean number of backlogged users Nb- 


A. Definition of Stability 

Let us study the conditions under which the described 
system is stable. In the followings we consider M and pq to be 
constant (i.e. stationary input). Once M and po defined, the 
expected channel input due to users in T state can be entirely 
described by the so called channel load line, which represents 
the relation between Gt and Nb for a certain communication 
scenario. For the finite population case, the channel load line 
can be defined as 


Gt 


M -Nb 


Po 


(4) 


while for M —oo the channel input can be described as a 
Poisson process with expected value A 0 so that Gt = A 
for any Nb, i.e. the expected channel input is constant and 
independent on the number of backlogged packets. 

Consider Figure representing various scenarios for 
CRDSA with Nf = 100 slots and Imax = 20. Equilibrium 
contours divide the (Nb,Gt) plane in two parts and each 
channel load line can have one or more intersections with 
the equilibrium contour. These intersections are referred to 


^Concerning th e va lues used for the Packet Loss Ratio, it is known from 
the literature jy (10) that the relation between PLR{Gjn) and G/jv can 
not be easily modelled in an analytical manner. For this reason PLR values 
used in this work are taken from simulations. 




(a) Stable channel (b) Unstable channel (finite M) 
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(c) Unstable channel (infinite M) (d) Overloaded channel 
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Fig. 3: Examples of stable and unstable channels for CRDSA with 
Nf = 100 and Imax = 20. Stable equilibrium points are marked 
with a black dot. 



























as equilibrium points. The rest of the points of the channel 
load line will belong to one of two sets: those on the 
left of the equilibrium contour represent points for which 
Gout > Gt, thus situations that yield to decrease of the 
backlogged population; those on the right represent points for 
which Gout < Gt, thus situations that yield to growth of 
the backlogged population. 

From the considerations above, we can gather that an 
intersection point where the channel load line enters the left 
part for increasing backlogged population corresponds to a 
stable equilibrium point, since it acts as a sink. In particular, 
if the intersection is the only one, the point is a globally stable 
equilibrium point (indicated as G^,N'g), while if more than 
one intersection is present, it is a locally stable equilibrium 
point (indicated as Gf ,iVj). If an intersection point enters the 
right part for increasing backlogged population, it is said to 
be an unstable equilibrium point (indicated as Gl^,Ng) in the 
sense that as soon as a statistical variation from the equilibrium 
point occurs, the communication will diverge in one of the 
two directions of the channel load with equal probability (as 
claimed in OJ ). 

Figure shows a stable channel. The globally stable 
equilibrium point can be referred as channel operating point 
in the sense that we expect the channel to operate around that 
point. With the word around we mean that due to statistical 
fluctuations, the actual Gt and Nb may differ from the 
expected value, however numerical results will show that 
averaging over the entire history of the transmission, values 


close to the expected ones are obtained. Figures 3b and 3c 


show unstable channels respectively for finite and infinite 
number of users. Analyzing this two figures for increasing 
number of backlogged packets, the first equilibrium point is a 
stable equilibrium point. Therefore the communication will 
tend to keep around it as for the stable equilibrium point 
in Figure 3 a and we can refer to it once again as channel 
operating point. However, this is not the only point of equi¬ 
librium since more intersections are present. Therefore, due 
to the abovementioned statistical fluctuations, the number of 
backlogged users could pass the second intersection and return 
to the right part of the plane, causing an unbounded increase 
of the expected number of backlogged users in the case of 
infinite M (Figure or an increase till a new intersection 
point is reached in the case of finite M (Figure [Tb|). In the latter 
case, this third intersection point is another stable equilibrium 
point known as channel saturation point, so called because it 
is a condition in which almost any user is in B state and 
Gout approaches zero. In the former case of infinite M, 
Nb will increase indefinitely and we can say that a channel 
saturation point is present for Nb oo. Notice that if M 
is finite, a stable channel can always be achieved using a 
sufficiently small value of Pr (the channel in Figure has 
the same parameters as the one in Figure 3b except for pr 
that is 0.5 instead of 1). However, when pr gets smaller, the 
corresponding average packet delay could get larger as we 
will show in Section V. Therefore a tradeoff between stability 
of the channel and average packet delay need to be faced 


in the design phase. Finally Figure shows the case of an 
overloaded channel. In this case there is only one equilibrium 
point corresponding to the channel saturation point. As for 
unstable channels with finite M, this channel can be rendered 
stable decreasing p^. Now, from the equilibrium contour we 
know that the point for the maximum expected throughput is 
^Qmax,Nmax'j, therefore the communication channel can be 
designed using Equation 0: 


G max 

T 


M - N^^^ 


Po 


(5) 


given that (M/Nf) ■ pq > G^ax must hold since the 
slope (known in the literature as m) of the channel load line 
must be a negative value. As an example, if po is a fixed 
and not modifiable design constraint, the number of users M 
that ensures maximum throughput can be calculated from 0. 
However the stability and the average packet delay need to be 
computed in order to verify they fulfill the design constraints. 
The former can be verihed as illustrated above, the latter can 
be computed as shown in the next section. 


IV. Packet Delay Model 

Assuming that a channel is operating at its channel op¬ 
erating point, we would like to know what is the expected 
distribution and mean delay associated to successfully trans¬ 
mitted packets. This can be entirely described using a discrete¬ 
time Markov chain with two states (Figure 0. The two 
states represent the state of a generic user that has a packet 
to transmit at the end of each frame. Therefore a frame 
duration is our discrete time unit for this Markov chain. The 
edges emanating from the states represent the state transi¬ 
tions occurring to users. These transitions depend on 
PLi?(GpjY, Nf,d, I max)- If a packet transmitting for the hrst 
time receives a positive acknowledgement (ACK) the user 
stays in T state while in case of negative acknowledgement 
(NACK) the user switches to B state until successful retrans¬ 
mission has been achieved and acknowledged. Therefore, here 
the packet delay Dpkt is considered as the number of frames 
that elapse from the beginning of the frame in which the packet 
was transmitted for the first time, till the end of the one in 
which the packet was correctly received. 



Fig. 4: Markov Chain for the Packet Delay analysis 

According to the definition given above, the delay distribu¬ 
tion is entirely described by 









Pr{Dpkt =n} = 


1 - PLR, 

PLR [pr (l-PLR)]- 
■[1 - Pr + PLR 


for n = 1 


for n > 1 


while the average expected delay is 
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0.65 


B 0.55 


Av[Dpkt] ='^n- Pr{Dpkt = n} 


Figure shows some examples of delay distribution. In this 
particular case, when pr gets smaller the distribution spreads 
over bigger values of delay and even though the probability 
that a packet is successfully transmitted at the first attempt 
increases, the expected average packet delay gets higher. This 
can be also explained noting that in this example, when 
decreasing pr the throughput decreases while the number of 
backlogged packets increases so that a packet will require 
more time to be successfully transmitted. Let us call 
the maximum throughput achievable. If M and po such 
that (M ■po)/Nf < Gq^^, the average packet delay at the 
channel operating point will always increase when decreasing 
Pr, since the equilibrium contour moves upwards. 
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Fig. 6: Simulated throughput for CRDSA with Nf = 100 slots, 
Imax = 20, Po = 0.143, Pr = 0.5, M = 350 


due to the aforementioned statistical variations, the throughput 
oscillates around a certain value. This value is the equilib¬ 
rium point. In fact the horizontal line represents the value 
Av[Gout\ = 0.486, that is the value obtained averaging the 
throughput over the entire simulation and is, as expected, a 
value really close to the one claimed in Figure Moreover, 
also the related average packet delay has been found to agree 
with the analytical results]^ 
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Fig. 5; Delay distribution for CRDSA with Nf = 100 slots, 
Imax = 20, M = 350 and po = 0.143 


V. Numerical Results 

In this section, the results of simulations are shown in order 
to validate the stability model described above. The simulator 
has been built according to the system description given 
in Section II, therefore perfect interference cancellation and 
channel estimation have been assumed. Moreover, neither the 
possibility of FEC nor power unbalance have been considered 
for our simulations. 

Figure shows the result of simulations for a communi¬ 
cation scenario with the same parameters as the example of 
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Fig. 7: Simulated throughput for CRDSA with Nf = 100 slots, 
Imax = 20, Po = 0.143, Pr = 1, M = 350 


Figure |7] illustrates an outcome of simulations for the 
communication scenario of the unstable channel described in 


stable channel described in Figure 3a It can be seen that 


Figure 3b In this case, the initial behaviour of the channel 
is such that the throughput oscillates around the channel 

^In the simulations only the average packet delay has been computed. In 
fact, regarding the delay distribution, a more complicated simulator is required 
in order to trace the identity of each packet so that each individual packet 
delay is known. Even though this can be accomplished, we have not built 
such a simulator yet. 


























































































operating point. However, differently from the previous exam¬ 
ple, this is not a globally stable equilibrium point. Therefore 
statistical fluctuations will sooner or later cause divergence 
from the channel operating point and the subsequent saturation 
of the channel with an average throughput that approaches 
zero. The time it takes for the channel to diverge from the 
channel operating point varies from simulation to simulation. 
In literature |j^ the expected time of divergence for Slotted 
Aloha is known as First Exit Time (FET). In our case this value 
is not easily computable since a Markov chain describing the 
communication and considering all the transition probabilities 
needs to be computed. This is computationally costly (as 
already claimed in GD) and we have not carried on methods 
and formulas for its calculation yet. 



Fig. 8; Simulated throughput for CRDSA with Nf = 100 slots, 
Imax = 20, A = 0.4, pr — 0.5, M —>■ oo lu the case before 
divergence from the operating point occurred 


Finally, Figure illustrates an outcome of simulations for 
the case illustrated in Figure 3c that is the case of infinite 
population. In particular, this outcome shows an occurrence 
in which divergence from the channel operating point has 
not occurred yet. As we can see, as long as the commu¬ 
nication takes place around the channel operating point, the 
same throughput and delay considerations as for the stable 
channel are valid. This example highlights that depending on 
the communication parameters, even though the channel is 
unstable, the FET could be so big that instability might be 
acceptable depending on application. However, if this is not the 
case, the communication will soon exit the stability region and 
the number of backlogged users will grow fast and indefinitely 
as shown in Figure 


VI. Conclusions 


In this paper a model for CRDSA based on equilibrium 
contour and extandable also to CRDSA-h- and IRSA has been 
proposed. This model allows to study and predict the stability 
and to calculate the expected throughput as well as the aver¬ 
age and expected distribution delay at the channel operating 
point, thus allowing a design that based on design constraints 
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Fig. 9: Number of backlogged users for CRDSA with Nf = 
100 slots, Imax = 20, A = 0.4, Pr — 0.5, M —>■ oo with divergence 
from the operating point after a certain time 


considers the tradeoff between stability, throughput and delay 
in order to provide an efficient and stable communication. As 
also suggested by reviewers (that we gratefully acknowledge 
for their suggestions and comments), future work will regard 
the analysis of a broader range of cases, the use of a different 
number of replicas than 2 in order to better understand this 
dependency and the study of such a communication scenario 
when introducing FEC codes as well as the possibility of 
packet power unbalance. 
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